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ABSTRACT 

We quantify the systematics in the size-luminosity relation of galaxies in the SDSS 
main sample (i.e. at z ~ 0.1) which arise from fitting different one- and two-component 
model profiles to the images. We use a novel argument to show that photometric infor- 
mation alone indicates that the majority of galaxies are not single-component systems. 
A SerExp model, a Sersic bulge plus exponential disk, for the two-components, is a 
better description than is a single Sersic profile; it is also better than the traditional 
deVaucouleurs bulge plus exponential disk model. In objects brighter than L*, fitting 
a single Sersic profile to what is really a two-component SerExp system leads to biases: 
the half-light radius is increasingly overestimated as n of the fitted single component 
increases; it is also overestimated at B/T~ 0.6. For such objects, the assumption of 
a single Sersic component is particularly misleading. However, the net effect on the 
size-luminosity relation is small, except for the most luminous tail. 

We then study how this relation depends on morphology. Our analysis is one of 
the first to use Bayesian-classifier derived weights, rather than hard cuts, to define 
morphology. Crudely, there appear to be only two relations: one for early-types, and 
the other for later-types (Es, SOs and Sa's are early-types, whereas Sbs and Scds are 
late). Closer inspection shows that within the early-type sample SOs tend to be 15% 
smaller than Es of the same luminosity, and, among faint late types, Sbs are more than 
25% smaller than Scds. Neither the early- nor the late-type relations are pure power- 
laws: both show significant curvature, which we quantify. This curvature confirms that 
two mass scales are special for early-type galaxies: = 3 x 10 10 M Q and 2 x 10 11 Mq. 
At M* > 2 x 10 11 M Q , the R — M* relation is the same as that of the bulge component. 
These same mass scales are also special for late types: there is almost no correlation 
between R and M* below the former, and almost no late-types above the latter. 

Our bulge-disk decompositions indicate that disks in both early- and late-types 
tend to be ~ 3x larger than bulges. Although the R — L and R — relations of 
disks of late types run parallel to those of the total light in late-types, these relations 
for the bulges of early-types show almost no curvature, and lie well below those for 
the total light in early- types especially at low L and M*. In fact, for narrow bins in 
velocity dispersion a, i?buigo oc £bulge (and similarly i?buigo oc M*buige): our SerExp 
bulges satisfy the virial theorem scaling, with bulges of fixed L but smaller a having 
larger sizes. The intrinsic scatter in the R — L relation decreases at large luminosity 
and/or stellar mass and should provide additional constraints on models of how the 
most massive galaxies formed. 

Key words: galaxies: structural parameters - galaxies: fundamental parameters - 
galaxies: evolution 



1 INTRODUCTION 
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The spatial (and color) distribution of star light in a galaxy 
is thought to encode information about its formation history. 
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In recent years, the correlation between size and luminosity 
for early-type galaxies has received much attention, because 
high redshift early-types appear to be more compact than 
their counterparts at low redshift (e.g. Trujillo et al. 2006; 
van Dokkum et al. 2008; Cimatti et al. 2008; Bruce et al. 
2012). However, both the size and the luminosity estimates 
are derived parameters, obtained by fitting to the observed 
surface brightness distribution. As a result, they depend on 
assumptions about the intrinsic shape of the surface bright- 
ness profile. E.g., if the fit assumes that galaxies are made 
up of two components or just one, and if two, whether they 
are modelled as the sum of an exponential and a deVau- 
couleurs (1948) profile, an exponential and a Sersic (1968), 
or two Sersics. 

The main goal of this paper is to quantify the system- 
atics on the local R—L relation which are associated with 
the choice of a particular model. In practice, 'local' means 
the 2 ~ 0.1 galaxies in DR7 of the SDSS Main Galaxy sam- 
ple ( Abazajian et al. 2009) . Because this sample is apparent 
magnitude limited, in practice, by R — L relation we always 
mean log 10 R fitted as a function of absolute magnitude (see 
Sheth & Bernardi 2012 for a simple description of the bias 
which would arise from fitting L as a function of R). And R 
denotes the radius which encloses half the total light L. (For 
exponential disks, this radius is 1.67 times the scale length 
of the exponential.) 

Our goal implies that we must fit the observed pro- 
files and determine the associated R — L relation using a 
variety of different models. We perform the fits using the 
PyMorph package, which can fit seeing convolved two compo- 
nents models to observed surface brightness profiles (Vikram 
et al. 2010). The algorithm is described and tested in Meert 
et al. (2012a, b). Tests on synthetic images show that when 
the fitted functional form is the same as the one used to 
generate the image, then PyMorph returns accurate values of 
the free-parameters (e.g., total light, half-light radius, bulge- 
total ratio). 

At first, we explore results from a single Sersic profile 
since the standard to date has been to use parameters from 
single Sersic fits. We will demonstrate that fits to a single 
Sersic profile are likely to return biased estimates of R and 
L, whereas fits to the sum of an exponential and Sersic pro- 
file should be less biased. Section [2] shows that the PyMorph 
derived R — L relation for single Sersic fits is offset from 
the standard in the literature (Shen et al. 2003). However, 
it is consistent with that obtained from a more recent de- 
termination of the Sersic parameters (Simard et al. 2011; 
hereafter Sll). Although we believe the Sll reductions are 
to be preferred to those in Shen et al., Appendix |A1 shows 
that they too are slightly biased - they appear to imply evo- 
lution in the n — L and R — L relations which we believe to 
be unphysical. This evolution is not present in the PyMorph 
reductions, so the remainder of the paper considers PyMorph 
exclusively. 

Section [3] compares the R — L relation based on single 
Sersic, Sersic + exponential, deVaucouleur + exponential, 
and single deVaucouleur fits, showing that the relations from 
single Sersic fits are offset to larger sizes and those from 
single deVaucouleur fits to smaller sizes, compared to those 
from the two-component fits. (In all cases, 'size' means the 
radius which contains half of the total light: in the case of 
two components, this is a complicated function of the light 



in each component and the scale radii.) Appendix IBI argues 
that although fits to a single Sersic profile are likely to return 
biased estimates of R and L, whereas fits to the sum of an 
exponential and Sersic profile should be less biased, the net 
effect on the derived R — L relation is small. On the other 
hand, fitting a realistic model is necessary to obtain sensible 
estimates of the intrinsic scatter around the mean R — L 
relation. This is the subject of Section [3.31 where we show 
that the intrinsic scatter correlates with a and decreases 
with increasing luminosity. 

Section 13.41 studies how the R — L relation depends 
on morphological type, using both the eye-ball classifica- 
tions of Fukugita et al. (2007; hereafter F07) and Nair et al. 
(2010; hereafter N10) and the Bayesian Automated Classifi- 
cations (hereafter BAC) of Huertas-Company et al. (2011). 
The latter are particularly interesting, because they are 
expressed as probabilistic weights - something we expect 
will become increasingly common in the next generation of 
large datasets. We explore the use of hard cuts based on 
these weights as indicators of morphology, as well as simply 
weighting each galaxy by the probability that it is one type 
or another. 

Bulge dominated galaxies do have disk components, and 
disk dominated galaxies have small bulges. In Section [4] we 
use our Sersic-exponential fits to study the R — L relation of 
the bulge components of the former and the disk components 
of the latter. In contrast to the curved relations we see for 
the total size-luminosity (size-M*) correlation, the R — L 
(R— M*) relation for bulges of early-types appears to be a 
power law, the slope of which is the same as that of the early- 
type R — L (R— M*) relation at fixed velocity dispersion. 
In addition, it shows that the two mass scales which are 
important for early-type galaxies, M* = 3 x 10 10 AfQ and 
2 x 1O 11 M0 (Bernardi et al. 2011a, b) are also special for 
late types. A final section summarizes our findings. 

When converting from angular sizes and apparent 
brightnesses to physical sizes and luminosities, we assume 
a flat ACDM model with £l m — 0.23 and a Hubble constant 
whose present value is Ho — 70 km s _1 Mpc _1 . 



2 THE SIZE-LUMINOSITY RELATION AT 
Z ~ 0.1 FROM SINGLE SERSIC FITS 

There is an analytic expression for the light enclosed within 
a given distance of the center of a circular Sersic profile. 
From this, the half light radius can be obtained easily. How- 
ever, if the object has axis ratio b/a 7^ 1, where b and a 
are the half-light radii along the principal axes of the im- 
age, then the corresponding expression must be integrated 
numerically. Since this can be time-consuming, it is usual 
to approximate this case by using the expression for a cir- 
cle, but with a suitably chosen effective circular radius. The 
most common choice is \/ba = a \Jbja, but Saglia et al. 
(2010) have recently shown that (b + a)/2 is more accurate: 
for bulge dominated systems the difference matters little, 
but it does matter for disks. Therefore, in the next subsec- 
tion (|2.1|l . where we compare with previous work, we use 
\fb~a. Thereafter, we use (b + a)/2. 
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Figure 1. The r-band singlc-Sersic based half-light radius (Rhl) versus total absolute magnitude (Mtot) relation for objects with n > 2.5 
(left) and n < 2.5 (right). In the panel on the left, our PyMorph determination is in good agreement with that based on single-Sersic 
parameters from Simard et al. (2011), but lies about 0.1 dex above, and is more curved than the fit reported by Shen et al. (2003). 
Symbols with error bars (joined by a solid curve for clarity) show the median half-light radius in bins of absolute magnitude. Dashed 
lines show the 16th and 84th percentile. In the panel on the right (objects with n < 2.5), except for the brightest objects, the Pymorph 
relation lies systematically 0.05 dex below that of Simard et al. 



2.1 Comparison with previous work 

Using the objects in an earlier SDSS data release, Shen et al. 
(2003) reported fits to the R — L relation for objects which 
had n > 2.5 and n < 2.5, where R, L and n were determined 
from a single Sersic fit to the light profile. Note that the 
Sersic parameters used by Shen et al. were estimated from a 
1-dimensional radial surface brightness profile (prof Mean), 
measured in ~ 5 — 10 azimuthally averaged annuli (Blan- 
ton et al. 2003). Thus, it is expected to be significantly less 
accurate than a 2-dimensional fit to the whole galaxy image. 

The Shen et al. relations for n > 2.5 and n < 2.5 are 
shown as the dashed and dotted lines in the left and right 
hand panels of Figure [T] respectively. The red and blue sym- 
bols with error bars show our determination of the single- 
Sersic based relation, where now R, L and n are from our 
PyMorph reductions, and the grey symbols and error bars 
show the R — L relation which follows from single-Sersic fits 
performed by Simard et al. (2011, hereafter Sll). 

For objects with n < 2.5, the Sll-derived relation runs 
parallel to that from Shen et al., but is offset to larger sizes 
by 0.05 dex, whereas the PyMorph-derived relation transi- 
tions from Shen et al. at low luminosities to Sll at high 
luminosities. For objects with n > 2.5 the PyMorph-derived 
relation lies about 0.1 dex above, and is more curved than 
the fit reported by Shen et al. The PyMorph and Sll based 
relations depart significantly from Shen et al. at the low 
and high luminous ends, where they curve upwards to larger 
sizes. For this reason, we are inclined to conclude that, at 
least at the bright end, Shen et al. is slightly biased. At the 
low end the curvature could be due to contamination by 
later-type galaxies. 

However, at the highest luminosities of objects with 
n > 2.5, the PyMorph and Sll relations are also slightly but 
significantly different from one another. Appendix IA1 shows 



that, in fact, at high luminosities, the derived magnitudes 
and sizes can be quite different: the correlated nature of 
these differences means that the R — L relation is only mod- 
erately affected. Appendix IA1 goes on to show that the Sll 
reductions appear to require rather dramatic evolution in n 
and R: both are larger at z = 0.2 than at z = 0.05. Since 
we believe this is unphysical, we conclude that the PyMorph 
reductions, which show no such systematic trend with z, are 
less biased, so we will use them in the remainder of this 
paper. 



3 THE SIZE-LUMINOSITY RELATION AT 
Z ~ 0.1 FROM ONE- AND 
TWO-COMPONENT FITS 

In this section we study how the R — L relation depends on 
the functional form for the surface brightness profile that 
was assumed when estimating R and L. We would espe- 
cially like to compare the effects of fitting one versus two- 
component models to the images. It is conventional to speak 
of these as 'bulge' and 'disk' components; while this is accu- 
rate for disk-dominated systems (typically later-type galax- 
ies), it may be better to think of the 'disk' component in 
bulge-dominated systems (typically early-type galaxies) as 
simply being a second component that is not necessarily a 
(thin, inclined) disk. 

We stated earlier that, for each component, we can ap- 
proximate the half-light radius by assuming the image is a 
circle of radius (6 + a)/2. But what should we do when we 
have two components? A natural choice would be to circu- 
larize each component using its own (6 + a)/2, and to then 
determine the half light radius of the sum of the circular- 
ized components, where each is weighted by the fraction of 



© 0000 RAS, MNRAS 000, 000-000 



4 Bernardi et al. 



Selection 


Ell 


SO 


Sa 


Sb 


Scd 


Selection 


Ell 


SO 


Sa 


Sb 


Scd 


EARLY- TYPES 












EARLY-TYPES 












Selected 












Selected 












P(E+S0)> 0.85 AND n > 3 
n > 2.5 


0.70 
0.44 


0.21 
0.18 


0.08 
0.20 


0.01 
0.13 




0.05 


P(E+S0)> 0.85 AND n > 3 
n > 2.5 


0.57 
0.32 


0.29 
0.23 


0.14 
0.29 



0.12 



0.03 


Missed 












Missed 












P(E+S0)< 0.85 OR n < 3 
n < 2.5 


0.10 
0.02 


0.43 
0.12 








P(E+S0)< 0.85 OR n < 3 
n < 2.5 


0.14 
0.02 


0.43 
0.07 








LATE- TYPES 












LATE- TYPES 












Selected 












Selected 












P(E+S0)< 0.15 AND n < 3 
n < 2.5 



0.01 


0.01 
0.04 


0.08 
0.11 


0.36 
0.33 


0.51 
0.45 


P(E+S0)< 0.15 AND n < 3 
n < 2.5 



0.01 



0.03 


0.15 
0.17 


0.41 
0.35 


0.39 
0.36 


Missed 












Missed 












P(E+S0)> 0.15 OR n > 3 
n > 2.5 






0.86 
0.75 


0.47 
0.41 


0.21 
0.16 


P(E+S0)> 0.15 OR n > 3 
n > 2.5 






0.79 
0.71 


0.35 
0.33 


0.17 
0.10 



Table 1. Eyeball morphological classifications from Fukugita et 
al. (2007). We set Ell (T = and 0.5), SO (T = 1), Sa (T = 1.5 
and 2), Sb (T = 2.5 and 3), and Scd (T = 3.5, 4, 4.5, 5, and 5.5). 

the total light that it contains (e.g. equation IC3|) . We have 
found that this approximation is quite accurate. 

3.1 Selection of morphological types based on 
hard-cuts 

Since the systematic biases may differ for bulge or disk dom- 
inated systems, we would like to separate out the effects of 
morphology from those on the functional form. Therefore, 
we will begin by first selecting objects of a single morpho- 
logical type. 

In practice this is difficult, because unambiguous deter- 
minations of the morphological type are not straightforward, 
although the task is slightly easier for bulge dominated sys- 
tems. We have chosen to select a sample of what we will call 
'early-types' on the basis of hard conservative cuts on two 
parameters which are available for each galaxy: the value 
of n obtained by fitting a single Sersic profile to the image, 
and the BAC probability p(E) that the object is an elliptical. 
We require n > 3 and p(E) > 0.85. These cuts by no means 
select all early-type galaxies; they are simply designed to se- 
lect a population which is very unlikely to be contaminated 
by later-types. Since our goal is to select objects of a single 
type, we are willing to sacrifice completeness for purity. 

In support of this assertion, Table Q] shows the mix 
of F07 morphological types in samples which are defined 
by hard cuts in the BAC p(type). Table [2] reports a simi- 
lar analysis which is based on the eye-ball classifications of 
N10 instead of F07. (Whereas BAC classifies galaxies into 
4 (E,S0,Sab,Scd) morphological types, N10 use the T-Type 
classification (—5 <T< 7) from the modified RC3 classifiers, 
and F07 use <T< 7 in steps of 0.5.) These Tables show 
that 91% and 86% of the resulting sample are indeed either 
Ell and SO. In contrast, requiring only n > 2.5 (as done in 
the past) yields a sample in which the E11+S0 fraction is just 
62% and 56% respectively. (A small fraction of the objects 



Table 2. Eyeball morphological classifications from Nair et al. 
(2010) who used T-Type classification using the modified RC3 
classifiers. We set Ell (T= -5 and T= -4), SO (T= -3, T= -2 
and T= -1), Sa (T= 0, T= 1 and T= 2), Sb (T= 3 and T= 4), 
and Scd (T= 5, T= 6 and T= 7). 

are Irregulars, which is why the numbers do not always add 
up to 100%.) Clearly our selection is much purer. As a mea- 
sure of its incompleteness, we also indicate the fraction of 
objects classified as Ells and SOs which do not make the cut. 
These fractions are 10% and 43% for the F07 classifications, 
and 14% and 43% for N10. This 'missed' fraction is much 
smaller if we only require n < 2.5, but we believe the price 
to pay in purity is unacceptable. 

The bottom halves of the two tables show a similar 
analysis of BAC and n cuts which are designed to produce 
a pure sample of later types. In this case requiring n < 3 
and p(E) < 0.15 yields a sample in which Sa + Sb + Scd 
account for 95% of the objects. If we only require n < 2.5, 
then Sa + Sb + Scd account for 89% of the objects so, for 
later types, the use of the BAC analysis does not make such 
a dramatic difference. 

3.2 Dependence of the R — L relation on model 
fitting and morphological type 

The panel on the left of Figure[2]shows the R~L relation ob- 
tained for this early-type sample (i.e. n > 3 and p(E)> 0.85) 
based on SDSS fits to a single deVaucouleurs profile; SDSS- 
based cmodel sizes defined by Bernardi et al. (2010), which 
are a crude combination of separate fits to a single expo- 
nential disk and a single deVaucouleurs profile; PyMorph fits 
to a two-component deVExp model; PyMorph fits to a two- 
component SerExp model; and PyMorph fits to a single Ser- 
sic profile. There are clear systematic differences between 
these relations, with the single Sersic and deVaucouleurs 
models returning the relations with the largest and smallest 
sizes, respectively. The various two-component based rela- 
tions are in good agreement except at the highest luminosi- 
ties (M r < —22), where the sample becomes increasingly 
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0.0242 
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ocrsH/xp ^(oa/oD/ocaj ) 
SersExp (P(Scd)) 


7.4437 
9.6010 
y.oloo 
7.8056 


0.8922 

1.0903 
1 m qo 

0.8396 


0.0266 
0.0311 
n hook 

0.0237 


SersExp (Early-type-Bulgcs) 
SersExp (Late-type-Disks) 


-2.0733 
6.4982 


0.0956 
0.6934 


0.0098 
0.0199 



Table 3. Luminosity-size relation. Early-types: p(Early-type)> 
0.85 and n > 3. Late-types: p(Early-type) < 0.15 and n < 3. 
Early-type-bulges: the bulge half-light radius versus the bulge lu- 
minosity for galaxies with p(Early-typc) > 0.85 and n > 3. Late- 
type-disks: the disk half-light radius versus the disk luminosity 
for galaxies with p(Early-type) < 0.15 and n < 3. 

contaminated by BCGs which are known to define steeper 
relations than the bulk of the population (e.g. Bernardi et 
al. 2007; Bernardi 2009). We argue in Appendix [B] that the 
SerExp reductions return less biased estimates of R and L. 

The panel on the right shows a similar analysis of the 
R — M* relation (stellar masses were computed as described 
in Bernardi et al. 2010 with a Chabrier IMF). Note that both 
R — L and R— M* are significantly curved. (Of course, if the 
stellar population models used to estimate M*/L are incor- 
rect, or if the IMF is mass-dependent, then this will modify 
the curvature in R — M*.) The curved dashed line shows a 
pure power law: the bottom panel shows that the deviation 
from this power law is substantial at \og 10 (M, /Mq) < 10.5 
and above log 10 (M* / Mq) > 11.3. These are the two mass 
scales identified by Bernardi et al. (2011a,b). 

We have repeated this analysis for a late-type sample, 
defined by requiring n < 3 and p(E) < 0.15. Although we 
do not show the corresponding plots here, we again see cur- 
vature. Rather, we illustrate this in Figure [3] which com- 
pares the SerExp-based R — L relation for our way of select- 
ing early- and late-type samples, with the more traditional 
cuts on n (larger or smaller than 2.5). The two ways of se- 
lecting the samples lead to very similar results, with the 
low luminosity early-types having smaller sizes, but defin- 
ing a steeper relation, so they would cross the R — L re- 
lation of late-types at about M r < —23 (beyond which 
there are few late-types anyway). We have also selected an 
intermediate-type population by requiring 2.5 < n < 3.5 
and 0.2 < p(E) < 0.4. Notice that this sample defines the 
same R — L relation as when we require our early-type se- 
lection (i.e., n > 3 and p(E) > 0.85), as well as that when 
we only require n > 2.5; we return to this in Section \3. 41 

We have quantified the curvature in these relations by 
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Figure 4. Top: Observed scatter around the mean (R\L) rela- 
tions for early-types based on fitting Sersic (solid red) and Ser- 
Exp models (solid green) to the images. Black solid curve shows 
the corresponding measurement from Shen et al. (2003). Dashed 
and dotted curves show a number of estimates from simulations 
of the measurement errors (see text for details). Grey dot-dashed 
line in top panel shows the rms difference between PyMorph and 
Sll sizes (both based on fitting a two-component SerExp). Bot- 
tom: Estimate of the intrinsic scatter around the Sersic (lower, 
red curve) and SerExp (upper, green) derived relations for early- 
types, obtained by subtracting in quadrature the red-dashed and 
green dotted curves from the corresponding red and green solid 
curves shown in the top panel. 

fitting to 

( lo giojJ^|°) =P0+PiO+ P2 O 2 ; (1) 

the coefficients of these fits for O — M r and O = M* are 
reported in Tables [3] and 0] 

3.3 Scatter in log(size) around the mean relation 

Our analysis allows us to make two interesting statements 
about the intrinsic scatter around the mean R — L relation 
for early-type galaxies. 

The top two jagged solid curves in the top panel of Fig- 
ure|4]show the measured scatter around the mean R—L rela- 
tion for SDSS early-types, when R and L are determined by 
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Figure 2. Dependence of derived size-luminosity (left panels) and size-stellar mass (right panels) correlations for early-type galaxies on 
the assumed surface brightness profile. Symbols with error bars (joined by a solid curve for clarity) show the median half-light radius 
in bins of absolute magnitude (left) and stellar mass (right). The SDSS fits to a single deVaucouleurs profile return a relation with the 
smallest sizes; our PyMorph fits to a single Sersic profile return the largest sizes. Of the relations which lie in between these two extremes, 
and which are almost indistinguishable at M < —21.5, the SDSS based cmodel sizes (defined by Bernardi et al. 2010) are the smallest; 
those based on our PyMorph fits to a two-component deVExp model are slightly larger; and those based on PyMorph fits to a SerExp 
model are largest. The curvature at the bright end appears to be due to an increasing incidence of BCGs, which define steeper relations 
(dotted lines) than the bulk of the early-type population. 
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Figure 3. Similar to previous figure, but now objects are selected using different hard cuts which define early-, late- or intermediate-type 
samples. Symbols with error bars (joined by a solid curve for clarity) show the median half-light radius in bins of absolute magnitude. 
Dashed lines show the 16th and 84th percentile. Note that this definition of intermediate's yields an R — L relation which is essentially 
the same as for the population with n > 2.5. 
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Sample/fit 


po 


Pi 


P2 


Sersic (Early-types) 
Sersic (Late-types) 
Sersic (n > 2.5) 
Sersic (n < 2.5) 


19.0933 
13.0054 
14.4995 
8.6098 


-3.9536 
-2.6438 
-3.1767 
-1.8301 


0.2070 
0.1393 
0.1742 
0.1014 


Sersic (P(E11)) 
Sersic (P(SO/Sa)) 
Sersic (P(Sa/Sb/Scd)) 
Sersic (P(Scd)) 


13.6593 
20.1092 
22.3082 
17.9815 


-2.9799 
-4.1549 
-4.4655 
-3.6102 


0.1635 
0.2166 
0.2275 
0.1862 


SerExp (Early-types) 
ScrExp (Late-types) 
SerExp (n > 2.5) 
SerExp (n < 2.5) 


13.4131 
11.2699 
12.5026 
9.5210 


-2.9324 
-2.3026 
-2.7875 
-1.9963 


0.1607 
0.1227 
0.1551 
0.1090 


ScrExp (P(E11)) 
SerExp (P(SO/Sa)) 
ocrrL/Xp ^r^oa/oD/ocaj ) 
SerExp (P(Scd)) 


12.8394 
19.2830 
lo.DloU 
11.7537 


-2.8246 
-3.9866 

-O. / 4z0 

-2.3957 


0.1557 
0.2079 

n i ooo 
u. lyzz 

0.1271 


SerExp (Early-type-bulges) 
SerExp (Late-type-disks) 


4.0853 
17.9763 


-1.4159 
-3.5683 


0.0992 
0.1831 



Table 4. Stellar mass-size relation. Early-types: p(Early-type)> 
0.85 and n > 3. Late-types: p(Early-type) < 0.15 and n < 3. 
Early-type-bulges: the bulge half-light radius versus the bulge 
stellar mass for galaxies with p(Early-type) > 0.85 and n > 3. 
Late-type-disks: the disk half-light radius versus the disk stellar 
mass for galaxies with p(Early-type) < 0.15 and n < 3. 

fits to a single Sersic (larger scatter) and to a SerExp model 
(lower scatter). This scatter is broader than the intrinsic 
one, because it includes a contribution from the measure- 
ment errors. For comparison, the smooth black curve shows 
the corresponding scatter reported by Shen et al. (2003). It 
is in reasonably good agreement with ours, except at the 
faint end, where we believe the enhanced scatter is due to 
increased contamination by spirals, for which the scatter is 
larger (as we show later). 

To estimate the intrinsic scatter, we must account for 
the broadening due to measurement errors. We estimate the 
errors on the sizes from fitting to the objects in the mock 
catalogs used in the Appendix B, where we know the input 
values. (See Meert et al. 2012a for details of how the mocks 
were generated.) The remaining dot-dashed curve shows the 
rms difference between PyMorph and Sll sizes returned by 
two-component SerExp fits to SDSS images, plotted as a 
function of the PyMorph SerExp absolute magnitude for the 
early-type sample. This is almost certainly an overestimate 
of the measurement error on the sizes; we have included it 
just to get a sense of the overall magnitude with which sys- 
tematic rather than random errors might affect the scatter 
in the R — L relation. 

The dotted and dashed lines show our simulation-based 
estimates of the measurement error on the sizes for an early- 
type sample. The lowest dotted line shows the rms scatter in 
log 10 R around the input value if the input profile is a single 
Sersic, and we fit it with a Sersic. In this, and all the cases 
which follow, we show this scatter as a function of the fitted 
(as opposed to the input) absolute magnitude. The other 
dotted line, which lies only slightly above the previous one, 
shows what happens if we fit a SerExp with a SerExp. These 
curves certainly underestimate the full measurement error, 
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Figure 5. At fixed velocity dispersion <r, the R — L relation is a 
power law whose slope is the same for all a, but whose zero-point 
increases as a decreases. 



since they are based on fits to smooth images, whereas real 
images may be lumpy, have spiral arms, etc. 

To get an idea of the magnitude of such effects, the 
two dashed curves show results from fitting a Sersic with 
a SerExp (lower) and a SerExp with a Sersic (upper). The 
differences between these and the dotted curves give an idea 
of the effect on the scatter of fitting an incorrect model to the 
data. The upper dashed curve is particularly interesting, in 
view of the fact that the SerExp model is more realistic (see 
Appendix B), whereas the single Sersic model is most often 
fit. Clearly, subtracting it in quadrature from the upper solid 
curve will lead to negative values at large luminosities. This 
is shown by the lower of the two curves in the bottom panel: 
at M r < —23 or so, the intrinsic scatter is consistent with 
zero. This, of course, does not mean that the R — L relation 
is intrinsically a line with negligible scatter. Rather, it is 
entirely a consequence of fitting an incorrect model. 

Recently, Nair et al. (2011) have used just such an argu- 
ment to claim that the R — L relation has no scatter. How- 
ever, their argument is based on Petrosian sizes and lumi- 
nosities; these are known to be inaccurate at large L, so the 
analysis above illustrates why their claim should be treated 
with skepticism. Indeed, the upper curve shows the result of 
subtracting (in quadrature) the upper dotted curve from the 
lower solid one, since both these are based on fitting to what 
we argued were more realistic models of the light profile. In 
this case, the intrinsic scatter is well-behaved: although it 
decreases steadily with M r , it does not go negative. 

Of course, since our estimate of the measurement error 
is really an underestimate, it is still possible that the intrin- 
sic scatter is smaller than we show. Therefore, we turn to 
what we believe is a much more effective way of showing 
that there is some intrinsic scatter. This method studies if 
the residuals from the relation correlate with other parame- 
ters, once correlations between the measurement errors have 
been accounted for. If they do, then there must be some in- 
trinsic scatter. 

Figure0shows the R—L relation for a number of narrow 
bins in velocity dispersion a. At fixed a, the R — L relation 
is a power law whose slope is 0.85 for all a but whose zero- 
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Figure 6. Size-luminosity (left) and size-M* (right) relations obtained by weighting objects by the BAC p(type). The results from 
SerExp fits are shown. The low L or M* part of the relation for Scds has the same slope as that reported by Shen et al. (2003) for 
their n < 2.5; and the intermediate L or M* part has the same slope they report for n > 2.5. Note that the relations for SOs are always 
indistinguishable from those for Es, and the Sab relations always lie between the E and Scd relations. Numbers in legend show the 
percentage of Ell, SO, Sa, Sb, Sbc and Irr galaxies classified by F07 with BAC P > 0.6. Using this selection we miss about 18% of Es, 
60% of SOs, 64% of Sab (37% Sa and 27% Sb) and 56% of Scd, respectively. 
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Figure 7. Comparison of the R — L relation in the morphologically defined samples of F07 (symbols and cyan curve) and N10 (magenta), 
with the fits defined by the BAC of Huertas-Company et al. (2011). All relations are in good agreement for E and SO galaxies; for 
comparison, the E relation is also shown in the other panels. F07 and N10 agree that Sa's define the same relation as Es and SOs, whereas 
Sb's are offset to larger sizes at smaller L. The HC-based results for Sab lie further from that for Es compared to those based on F07 
and N10 for Sa's, but are in good agreement for Sb's; however, for Scd's they lie closer to the E relation than do F07 or N10. 
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point increases as a decreases. In effect, this shows clearly 
that the scatter around the mean R — L relation correlates 
with cr; it is not all due to measurement errors. The slope of 
0.85 is consistent with previous work (Bernardi et al. 2003; 
Bernardi 2009); while steeper than the slope of 0.64 associ- 
ated with averaging over all a, it is less than unity - a fact 
we return to later, when we discuss the bulges of early-types 
(see Figure [TO)) . 

We end this subsection with the observation that the 
intrinsic scatter appears to be smallest for the most lumi- 
nous objects. Since it is commonly believed that mergers will 
affect the scatter of scaling relations such as this one, our 
overestimate of the intrisic scatter in the R — L relation pro- 
vides a new constraint on models of how the most massive 
galaxies must have formed. E.g., Shen et al. (2003) argue 
that many minor mergers may be more consistent with the 
shape and scatter of the R — L relation than are few ma- 
jor mergers. Other work has also explored constraints which 
come from the scatter (Shankar et al. 2012); it will be inter- 
esting to revisit this question in light of the mass-dependence 
we believe we see. 



3.4 Dependence on morphology based on the 
Bayesian Automated Classifier 

In the previous section we used a hard cut on the BAC 
probability to determine morphology. Since this is not quite 
in the spirit of why such probabilities were generated in 
the first place, Figure [6] shows the size-luminosity (left) and 
size-M, (right) relations obtained by weighting objects by 
p(type) as determined by BAC. The best-fit parameters of 
equation flj to these relations are reported in Table [3] and 

El 

Notice that the Sab's define a relation which lies be- 
tween that defined by Scd's on the one hand and Ell and 
SOs on the other. However, the Sab class is difficult to de- 
fine (c.f. Tables [JJ and f2|. Figure [S] reports the percentage of 
Ell, SO, Sa, Sb, Sbc and Irr galaxies classified by F07 with 
the BAC p > 0.6 - note that this is a different cut from 
that used in Table [JJ More than a third of the objects with 
p(S0) > 0.6 are Sa's, and about a fifth of the objects with 
p(Scd) > 0.6 are Sb's. Conversely, of the objects which have 
p(Sab) > 0.6, about one third are Scd's. 

To address this more closely, Figure [7| shows the R — L 
relations in the F07 eye-ball classified subsamples. The cyan 
curves show fits to these subsamples, and the magenta curves 
show fits based on the N10 (eyeball) classifications. The two 
are in quite good agreement. To emphasize the fact that 
the relation is different for the different subsamples, the red 
solid curve, which is the same in each panel, shows the BAC- 
based relation for p(Ell). The orange, green and blue curves 
(in the relevant panels) show the BAC-based relations for 
p(S0), p(Sab) and p(Scd). These are in good agreement with 
the F07 and N10 based relations for E and SO galaxies. 

Note that F07 and N10 agree that Sa's define the same 
relation as Es and SOs, whereas Sb's are offset to larger sizes 
at smaller L. This suggests that combining Sa's and Sb's into 
a single type may be problematic. Indeed, the BAC-based 
results for Sab lie further from that for Es compared to those 
based on F07 and N10 for Sa's, but are in good agreement 
for Sb's; however, for Scd's they lie closer to the E relation 
than do F07 or N10. These small but systematic differences 




Figure 8. At fixed luminosity, Es tend to be about 0.06 dex 
larger than SOs, although this offset depends slightly on how R 
and L were determined. 



between the BAC and eye-ball based results suggest that 
combining Sab's into a single class results in a weighted sum 
of the relations defined by E's and Scd's. 

Figure [7] shows that the curvature in the R — L relation 
is such that, for Scds, there is almost no correlation at Mr > 
—20.5. This flattening at low luminosities is also evident for 
the other morphological types, and is more pronounced in 
the R—M* relation shown in the right hand panel of Figure[6] 
(see also Figure [9] below) . Indeed, Figure [6] shows that at 
log 10 M*/Mq < 10.5, even the samples weighted by p(Ell) 
and p(S0) tend to have essentially no correlation between R 
and M, . 

This is the same mass scale at which a number of 
other early-type galaxy scaling relations change qualitatively 
(Bernardi et al. 2011a,b). Since Figures [6] and [7] indicate that 
it also appears to be significant for late-type galaxies, it is in- 
teresting to ask if the other mass scale identified by Bernardi 
et al., M* = 2 x 1O 11 M0, is also significant for late-types. 
Figure [6] shows that, in fact, this mass scale seems to set the 
limit above which there are essentially no late-type galaxies. 
Figure [7] tells a consistent story: although there are many Es 
and SOs brighter than M r — —23, there are no Sa, Sb or Scds 
with luminosities this large. Bernardi et al. suggested that 
this mass scale was associated with merger histories that 
were dominated by major dry mergers; since such mergers 
would destroy disks, the fact that we see no late-types above 
this mass scale is, perhaps, not surprising. 



3.5 Small but statistically significant difference 
between Ellipticals and SOs 

Above, we noted that there is essentially one R — L rela- 
tion for E, SO and Sa galaxies. However, our sample is large 
enough to detect small but significant differences within the 
early-type (E and SO) sample. A closer look at Figures [6] 
and [7] indicates that SOs are slightly smaller than Es of the 
same luminosity. Figure [8] shows that this offset is about 
0.06 dex, although it depends slightly on how R and L were 
determined. This is particularly interesting in view of recent 
work at z ~ 1, based on the Sll reductions, which shows a 
similar offset of about 15% for the SDSS sample growing to 
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Figure 9. Similar to Figure [3] but now contrasting the R — L 
and R — M* relations for early- types with that for their bulges, 
and the relation for late-types with that for their disks. 

~ 40% at z ~ 1 (Huertas-Company et al. 2012). Both the 
sign of the trend and its evolution deserve further study. 



4 BULGES AND DISKS 

One of the virtues of our SerExp decompositions is that 
it allows us to study the scaling laws of disks and bulges. 
Figure [§] contrasts the R — L relation for early-types with 
that for their bulges, and the relation for late-types with 
that for their disks. 

The i?disk — idisk relation runs parallel to the R — L 
relation for late-types; 7?disk tends to be 0.1 dex larger than 
Rhi- That iidisk > Rhi is not surprising, since we know that 
late-type galaxies host small bulges which will contribute to 
the light at small radii. But that this should have produced 
a constant offset is not obvious. We address this question 
shortly. 

The bulges are more interesting. In contrast to when 
the total light was used, there is almost no curvature in the 
relation for bulges. It is well approximated by a single power- 
law: (.Rbuigcl J^buigo) oc -t/buigc- The amplitude of the power 
law is such that the relation for bulges is approximately 
the same as for the total at very large luminosities; as L 



Figure 10. Same as Figure [5] but showing iJbulgc ~ ^bulge for 
a number of bins in total velocity dispersion a. Replacing Lbulgo 
with M, of the bulge yields the same result: The relation is a 
power law whose slope is 1 for all a, but whose zero-point increases 
as a decreases. 

decreases, the R — L relation curves away from the i?bui gc — 
ibuige relation, towards larger sizes. 

The power-law nature of the bulge relation suggests a 
picture in which the curvature in the early-type R — L rela- 
tion arises as a consequence of adding a disk component to 
bulges. However, there is an interesting puzzle: recall that 
Figure \S\ shows the R — L relation for a few narrow bins 
in velocity dispersion. This relation also has no curvature; 
remarkably, it runs parallel to the 7?bui gc — ibuigc relation, 
having slope ~ 0.85. To explore this further, Figure fTUl shows 
the analogue of Figure [5] the J?bui gc — ibuigc for fixed bins in 
a. In this case, the slope is 1. Replacing Lbuigo with M*buige 
makes no difference. I.e., our SerExp bulges exhibit the scal- 
ing expected from the virial theorem. 

The bottom panel of Figure [9] shows the corresponding 
relations when plotted as a function of A/*. The J?disk — 
Af*disk relation is again offset from that for all late-types; 
both are slightly more curved than their counterparts in the 
panel above. At log 10 M*/Mo < 10.5, the flattening of the 
relation with respect to the slope at large M* is such that 
there is almost no correlation between J?disk and M*disk- 
This flatness at the faint, low mass end is similar to that for 
Scds (see Figures [6] and [7] and related discussion). 

The bottom panel also shows that the J?bui g c — Af*buige 
relation sits on top of that for early types at the largest 
masses, suggesting that the second component which con- 
tributes somewhat to the light contributes little to the mass. 
It is worth noting that this happens at the same mass scale, 
M» = 2 x 10 n Af Q , which Bernardi et al. (2011a) noted was 
significant for early-types, and above which there appear to 
be no late-type galaxies (as is clear from this figure, as well 
as from Figures [6] and [7)l . 

Of course, as we cautioned before, the conversion from 
L to M* depends on M*/L, which in turn depends on stellar 
population modelling as well as on an assumption about how 
the IMF depends on galaxy mass, so these assumptions will 
affect how this measured curvature - and this mass scale in 
particular - should be interpretted. There is another reason 
to be cautious: our M* estimates assume that M*/L for 
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Figure 11. Correlation between -Rhi/-Rbulgc and B/T for early- 
types (top panel) and between i?hl/^?disk an d B/T for late- types 
(bottom panel). Although only objects with —21.5 > M r > —22.5 
are shown, we see qualitatively similar behaviour at other lumi- 
nosities. Dashed, dotted and solid curves show the expected scal- 
ing for n = 4 bulges with exponential disks having i?disk/^bulge = 
2,4 and 6. 

the individual components is the same as that for the total. 
Since we are looking at bulges of early types and disks of 
late types this assumption, while crude, should not be wildly 
wrong (the same would not be true for, e.g., the bulges of late 
types). Nevertheless, one might imagine that, as a result, 
we slightly under(over)-estimate the mass in the bulge(disk) 
component. 

Why are these relations for the bulge and disk compo- 
nents so different from those for the total light? 

4.1 The smallness of bulges 

Suppose we start from the power-law iibulge - ibulge relation, 
ibuige, with a given value of B/T, specifies a total magni- 
tude Mbuigc + 2.5 log 10 (B/T), for which the associated half- 
light radius would be given by log 10 i?bui g e+log 10 Rh/Rbuigc- 
The dots in the top panel in Figure [TT] show the shift in 
size which is associated with each shift in magnitude (it- 
self due to B/T), for SerExp fits to the early- type sample 
with —21.5 > M r > —22.5. There is clearly a strong correla- 
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Figure 12. Dependence of the ratio of disk to bulge size on total 
luminosity for early- (red) and late-type (blue) galaxies. 



tion between i?hi/-Rbui g c and B/T (at this fixed M r ) for the 
early-type sample. If the total is 0.55 mags brighter than the 
bulge (B/T= 0.6), then the half-light radius of the total is 
about 0.35 dex larger than that of the bulge. What causes 
this? 

The curves show the expected relations for a deVau- 
couleur bulge with exponential disk. These depend on the 
ratio -Rdisk/-Rbui go (we show 2, 4 and 6) but they are indepen- 
dent of the total luminosity (Appendix [C] shows why). The 
dependence on n is weak (repeating the analysis with n — 3 
brings the curves into better agreement with the measure- 
ments; n = 6 shifts in the opposite direction). Matching the 
data indicates that -Rdisk/flWgc ~ 5 at B/T < 0.7, suggest- 
ing that the correlation is caused by the fact that PyMorph 
uses disks with rather large scale lengths to account for the 
fact that a Sersic bulge is not, by itself, always a good match. 

Whether or not these large scale lengths are physically 
reasonable is an open question, but we show in Appendix [Cl 
that these tend to be objects for which the single Sersic 
fit returns large values of n > 5; these extended second 
components do appear to be necessary to provide a good 
fit. Indeed, fitting SerExp images with a single Sersic profile 
requires large values of n if 0.4 <B/T< 0.7 (bottom left 
panel of Figure IB1[) . 

The bottom panel shows a similar analysis of the late- 
type sample: -Rhi/-Rdisk as a function of (1-B/T). Most of 
the sample has B/T < 0.2 for which log 10 (i?hi/-Rdisk) dif- 
fers from zero by —0.05 dex or less. Although this is in the 
opposite direction to the shifts for early-types (as it should 
be), the resulting estimate of -Rdisk/-Rbui gc ~ 5 is similar. Of 
course, in this case, we expect J?disk 3> -Rbuigc, so the value 
of 5 does not require further explanation. 

We can, of course, directly measure the ratio 
^disk/^buigo for the objects in our early- and late-type sam- 
ples. Figure [12] shows that this ratio is indeed large, with 
only a weak dependence on L, and a somewhat larger scat- 
ter for early-types. The actual median value, ~ 3 — 4, is 
slightly smaller than the value of 5 we derived from the pre- 
vious figure on the basis of the idealization that all galaxies 
were deVaucouleur bulges with exponential disks. Hence, we 
conclude that the differences between the relations shown in 
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Figure 13. Top: Observed (top) and intrinsic scatter (bottom) 
around various R — L relations as labelled (format similar to 
Figure [4} . In all cases, our upper limit to the intrinsic scatter 
decreases at large luminosities; this is particularly dramatic for 
later- type galaxies. 



Figure [9] can be traced to the fact that bulges are substan- 
tially smaller than disks. 



4.2 Scatter 

Before ending this section, Figure [F3l shows our estimate of 
the measured and intrinsic scatter around the mean R — L 
relations defined by bulges and disks, and compares them 
with corresponding estimates for early-types and late-types. 
Notice that the measured scatter is substantially smaller 
around the early-type relation than around any of the oth- 
ers. Since we argued earlier that the Shen et al. (2003) early- 
type sample is contaminated by later-types, we believe this 
explains the difference between their results and ours in Fig- 
ure [4] Note also that the scatter around the relation for 
bulges is substantially larger than for the others. 

Our estimates of the intrinsic scatter (shown in the 
bottom panel) come from subtracting, in quadrature, the 
measurement errors seen in simulations (dashed lines) from 
the total scatter measured in the data (corresponding solid 
lines), following the method described in Section EOl For this 
reason, we are almost certainly overestimating the intrin- 



sic scatter. Nevertheless, it is interesting that for late-types, 
disks and bulges, our estimates indicate that the intrinsic 
scatter decreases at large luminosities. For early-types this 
decrease is less dramatic, with the scatter perhaps even lev- 
elling out at large luminosities. We believe these differences, 
along with the power-law nature of the bulge R — L relation, 
will prove to be useful for improving our understanding of 
how massive galaxies assembled their mass (e.g. Shankar & 
Bernardi 2009; Shankar et al. 2010). 



5 SUMMARY 

We used our automated image decomposition algorithm 
PyMorph to study the effects of systematics in the size- 
luminosity relation of galaxies in the SDSS main sample (i.e. 
at z ~ 0.1) which arise from fitting different models to the 
images. 

In Appendix A we argued that PyMorph returns more 
physically reasonable results than does the algorithm of Sll 
(e.g. Figures IA2I and IA3I and related discussion). And in 
Appendix B we showed that SDSS photometric information 
alone indicates that the majority of galaxies are not single- 
component systems, but have (at least) two-components. 
These are better modeled as a Sersic bulge plus exponen- 
tial disk, rather than the traditional deVaucouleurs bulge 
plus exponential disk. 

For objects brighter than L,, the commonly adopted 
procedure - of fitting a single Sersic profile to what is really 
a two-component SerExp system - leads to biases. The half- 
light radius is increasingly overestimated as n of the fitted 
single component increases; it is also overestimated around 
B/TscrExp ~ 0.6. For such objects, the assumption of a sin- 
gle Sersic component is particularly bad. However, the net 
effect on the size-luminosity relation is small, except for the 
most luminous tail (Figure [2Jl. 

On the other hand, fitting a realistic model is necessary 
to obtain sensible estimates of the intrinsic scatter around 
the mean R—L relation. Having done this, we showed that 
the scatter in sizes correlates with velocity dispersion, and 
the rms scatter decreases at large luminosity (Figure [4]), al- 
though for early-types it may level off to a constant value 
of about 0.1 dex at large luminosities. This should provide 
tight constraints on the nature and number of mergers re- 
quired to assemble the most massive galaxies. 

Our Figure [S] shows one of the first use of Bayesian 
classifier-based weights in the estimation of the R — L scaling 
relation for different morphologies (e.g. Aguerri et al. 2012). 
We found that, even if we allow for finer bins in morphology, 
there seem to be only two fundamental R — L relations, both 
of which are slightly but statistically significantly curved 
(Figures [3] [6] and [7] and Tables [3] and [4} . Of course, a closer 
inspection does reveal subtle dependences on morphology. 
Amongst early-types, SOs tend to be about 0.06 dex smaller 
than Es of the same luminosity (Figure [Sj . This difference 
smaller than the ~ 40% reported by Huertas-Company et 
al. (2012) at z ~ 1. This subtle difference between Es and 
SOs is particularly interesting in view of the fact that the 
two types show very different trends as a function of age 
(Bernardi et al. 2010), so we expect that it, and its evolu- 
tion, should yield interesting new constraints on models of 
how early-type galaxies assembled their stellar mass. Sim- 
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ilarly, amongst late-types, faint Sbs tend to be ~ 0.1 dex 
smaller than Scds of the same luminosity, but these differ- 
ences decrease as luminosity increases. 

Our two-component fits allowed us to study the R — L 
relations for the bulge and disk components themselves. Al- 
though the R — L relations for the total light in early- and 
late-types are curved, the relation defined by the bulges in 
early-types is remarkably straight: (-Rbui g e|£buigc) oc Lbulge 
(Figure [9]). The relation for disks runs parallel to the R — L 
relation for late type galaxies, being offset upwards by about 
0.1 dex. For disks, this curvature is so pronounced that, at 
the faint, low mass end, there is almost no correlation be- 
tween R and L or M* (Figures [Jj and [9]). We argued that, 
both for early and late type galaxies, these differences arise 
because PyMorph uses disk-components for which the half 
light radius is ~ 3 — 4 times larger than that of the bulge 
(Figures 1 1111 121 and Appendix [Cl> . It is not clear if for early- 
types this is physically reasonable - but extended second 
components are clearly necessary for the SerExp fits (Fig- 
ures [C2] and [C3]) . 

The two mass scales, M* = 3 X 1O 1O M and M» = 2 x 
lO n M0, previously identified by Bernardi et al. (2011a,b), 
also appear in Figure [5] For early-types, the former, is, 
among other things, the mass scale at which galaxies are 
maximally dense. Below this scale the R—M* relation curves 
upwards with respect to the power law which best describes 
the full range of M* (Figures [2] and [6]). Bernardi et al. sug- 
gest that this is because the disk component becomes more 
significant at these low masses. 

The larger mass scale is where the R — L relation of 
early-types curves upwards with respect to the power law 
which best describes the full range of M*. Bernardi et al. 
attribute this to a change in the assembly histories - to 
ones in which major dry mergers become important. So it is 
interesting that we find that it is at this mass scale that the 
bulge and total R — M* relations become the same, despite 
being very different at smaller masses (Figure [9]). This is 
particularly remarkable in light of recent work showing that 
early-types below this mass scale tend to be fast rotators 
(Cappellari et al. 2012). It may be that our SerExp bulge- 
disk decompositions of the images are reflecting this change 
in the kinematics. 

Our analysis indicates that these same two mass scales 
are also significant for late-type galaxies. At M* < 3 x 
10 10 Mq, the R — M* relation for late-types (and their 
disks) flattens significantly (Figures [5J [7] and O ; and M* = 
2 x 10 Mq marks the mass scale above which there are 
almost no late-types (Figures [Jj and [9]). 

Given the large differences between the relation for 
bulges and that for early-types at smaller masses and lu- 
minosities (Figure [9]), it is remarkable that the slope of the 
R — L relation for bulges is essentially the same as that for 
early-types within a fixed bin in velocity dispersion (Fig- 
ure [5]). Even more remarkable is the fact that, at fixed a, 
ibuigc oc i?buigc (Figure [TO)) : i.e., our SerExp bulges ex- 
hibit the scaling expected from the virial theorem. Why this 
should be so is an open question, but it does suggest that de- 
partures from the virial theorem scalings for the total light 
are entirely due to the presence of a disk component. 

Finally, we find that the scatter around the mean R — 
L relation decreases as L increases (and similarly for R — 
M*), except for early-types, where it may flatten at 0.1 dex 



(Figure I13[) . We expect this to provide a useful probe of 
how massive galaxies assembled their mass (e.g. Shankar & 
Bernardi 2009; Shankar et al. 2010). 
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Figure Al. Differences between the Sll and PyMorph reductions tend to be of order 0.04 mags fainter or 0.03 dex smaller in size, 
except for M r < —22.5 where PyMorph tends to be bigger and brighter if n > 2.5. 



APPENDIX A: SYSTEMATIC EFFECTS IN 
THE SIMARD ET AL REDUCTIONS 



The main text showed that the R — L relation from single- 
Sersic fits using PyMorph is in reasonably good agreement 
with that based on parameters from Simard et al. (2011). 
However, Figure lATl shows that, although the two algorithms 
return similar sizes and luminosities for objects with n < 2.5 
(PyMorph is about 0.03 dex smaller and 0.03 mags fainter), 
the PyMorph sizes and luminosities are systematically larger 
at large M to t- This bias for the biggest galaxies is particu- 
larly evident when shown as a function of Mp y Morph- 

Since the R— L relation of the largest galaxies is particu- 
larly timely, we would like to determine which reductions are 
more reliable. Figures [A2I and [ A3I show that the Sll reduc- 
tions indicate substantial recent evolution toward smaller 
n and R at fixed L especially at larger L. We believe this 
evolution is unphysical, so conclude that the Sll reductions 
suffer from systematic biases. No such evolution is seen in 
the PyMorph reductions, so we use them exclusively in the 
main text. 



APPENDIX B: SERSIC INDEX AND B/T 
RATIO IN SDSS GALAXIES: EVIDENCE FOR 
TWO COMPONENTS IN THE SURFACE 
BRIGHTNESS PROFILE 

The main text considered the size-luminosity relation, not- 
ing that the derived relation depends systematically on the 
assumed form of the surface brightness profile. In this Ap- 
pendix, we provide an analysis of the light profiles of SDSS 
galaxies which we believe strongly suggests that fitting to a 
SerExp model returns the least biased answers. 



Bl How many components? 

As noted in the introduction, there has been considerable 
interest in developing accurate descriptions of the projected 
surface brightness distribution of galaxies. 

One approach to this problem is to fit the free param- 
eters of a predetermined functional form to the observed 
surface brightness profile. These derived free parameters 
(typically, these are expressed in terms of the scale which 
contains half the total light, and the surface brightness at 
this scale) are more useful if the functional form itself actu- 
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Figure A2. Our determination of the n — L relation (symbols connected by solid lines) shows little or no rcdshift dependence. The 
sudden drop in n at the faint end of each redshift sample is due to the bimodal distribution in n at each L; it has nothing to do with 
evolution. Except for this, our determination shows little or no redshift dependence; in contrast, at high luminosities, the Simard et 
al. reductions appear to shift systematically towards smaller n as redshift decreases. We believe the implied evolution is unphysical, so 
conclude that the Simard et al. reductions are systematically biased. 



ally does provide a good fit to the profile. A simple version 
of this approach is to fit many different functional forms 
to the data, and then select the one which provides the 
best fit (in some suitably quantified way). For example, the 
Sloan Digital Sky Survey (SDSS; Stoughton et al. 2002) fits 
both exponential (1(0) oc exp(— 0/0i)) and de Vaucouleurs 
(oc exp[— (d/dt) 1 ^ 4 }) profiles to the image, along with an es- 
timate of which fits better. 

The exponential and de Vaucouleurs (1948) profiles are 
special cases (n = 1 and 4) of the Sersic (1968) profile 
(oc exp[—(9/9 n ) 1 ^"]). With sufficiently good data, it is pos- 
sible to simply fit a Sersic profile to the data, leaving the 
fitting procedure to determine n as well. If galaxies really 
are intrinsically single Sersics with a wide range of n, then 
the parameters (e.g. half-light radius) returned by forcing 
n — 1 or 4 in the single component fits will generally be bi- 
ased. Across the population as a whole, the derived value of 
n spans a wide range, sometimes being as large as ~ 8 or 10 
(e.g. Simard et al. 2011 and references therein), suggesting 
that forcing n — 1 or 4 is ill-advised. 



Of course, it is not obvious that the light profile should 
be fit using a single component. The stellar kinematics in 
many galaxies indicate that the stars define more than one 
dynamical component. Examples include counter-rotating 
disks, as well as disk systems with bulges or bars in their 
centers (e.g. Emsellem et al. 2011). Evidence for more than 
one component is often seen in the chemical composition as 
well (e.g. Johnston et al. 2012). In such galaxies, it is inter- 
esting to see if the light profile also indicates the presence 
of more than one component. 

This has motivated studies which model the observed 
profile as the sum of an exponential and a de Vaucouleurs 
profile; what we will call the deVExp model. (Of course, since 
there are now more free parameters to be fit, better, higher 
resolution data are required. In this context, it is worth not- 
ing that Sersic's initial motivation was to fit a functional 
form with fewer free parameters which would allow one to 
interpolate between two-component systems having varying 
fractions of an n = 4 bulge and an n — 1 disk.) It is common 
to report the result of such two-component fits in terms of 
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Figure A3. Similar to Figure lA2l but now for the R — L relation: little or no redshift dependence is seen in our sample; in contrast, 
at high luminosities, the Simard et al. reductions (symbols connected by dashed lines) imply evolution towards smaller sizes as redshift 
decreases. We believe this implied evolution is unphysical, so conclude that the Simard et al. reductions suffer from systematic biases. 



the fraction of the total light that is in the bulge (de Vau- 
couleurs) component: B/T. Correlations of these B/T val- 
ues with other parameters (e.g. luminosity) are then used to 
constrain formation history scenarios. 

On the other hand, if galaxies really are single compo- 
nent Sersics, and one attempts to fit them with two compo- 
nent deVExp profiles, then one will infer an entirely spurious 
B/T value (the profile was, after all, just a single compo- 
nent). This spurious B/T will correlate with other parame- 
ters if n itself does, complicating the interpretation of such 
correlations. Indeed, some have argued that the evidence for 
two-components in the light profile is sometimes just a con- 
sequence of trying to fit what is really a single component 
Sersic with a linear combination of exponential and deVau- 
couleurs profiles. (This leaves unanswered the question of 
why dynamically or chemically distinct components do not 
leave a signature in the light.) 

One way to address this question is to fit the image 
with the sum of two Sersic profiles, each with its own value 
of n, and then see if allowing for the second component does 
indeed provide a statistically significant improvement in the 
accuracy of the fit (once one accounts for the increase in 



the number of fitted parameters). In what follows we will 
perform a slightly simpler version of this: we force one of 
the components to have n = 1, while leaving the other to 
be determined by the fitting procedure. We then provide 
a novel argument which indicates that this SerExp model 
is indeed a better approximation to the surface brightness 
profiles of real galaxies than is either a single Sersic, or the 
deVExp model. 

We are not the first to have come to this conclusion; 
e.g. Allen et al. (2006) argued that at least half of the ~ 10 4 
galaxies at z ~ 0.1 in the Millenium Galaxy Catalog are two 
component SerExp systems, and Simard et al. (2011) have 
recently performed a similar analysis of ~ 10 6 SDSS galax- 
ies. But our argument for why we believe two components 
are needed is new. 

To gain intuition, section [B2l shows the result of fitting a 
variety of synthetic images (generated using either a single 
or two-component models) with single Sersic, deVExp and 
SerExp profiles. This section also presents a similar analysis 
of real galaxies selected from the SDSS. Section lB3l discusses 
some biases which arise from fitting the image with a single 
Sersic. 
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Figure Bl. Fitted n ser vs fitted B/T for simulated images which were generated using a single component Sersic profiles (top), or 
two-component deVExp (middle) or SerExp profiles (bottom). The two left columns show n ser , returned by fitting a single Sersic profile 
to the image, versus B/T, returned from fitting a deVExp profile; the two right columns show the same n sor , but now B/T comes from a 
SerExp fit. For each pair of columns the left column shows the density across the full sample, while the right column shows the density 
for four bins, colored by input n ser (top; red showing larger n BCT ) and input B/T (middle and bottom; red showing larger B/T). 




Figure B2. Similar to Fig. IB 1 1 but for real galaxies. Fitted (single component) n ser vs fitted B/T using the two component deVExp 
fit (two left columns) and the two component SerExp fit (two right columns). Colors represent the probability that the galaxy is an 
early- type (red is highest probability). 



B2 Fits to synthetic images 

In this section we show the result of using PyMorph (Meert 
et al. 2012a,b) to fit a variety of synthetic (mock) and real 
(SDSS) galaxies. We contrast what happens when PyMorph 
is forced to fit an image using only a single Sersic compo- 
nent, to when it is allowed to use two Sersic components, 
one with n = 1 and the other free: the SerExp model. For 
the two-component fits, we first show results when n of the 
Sersic component is set to 4, since this corresponds to the 
traditional 'deVaucouleurs bulge + exponential disk' deVExp 
fits, and then when n is allowed to be a free parameter, de- 
termined by the fit. 

In all the results which follow, the parent distribution is 
essentially a random subset of the SDSS DR7 main galaxy 
sample, which is magnitude limited to m r < 17.7. We fit 
each object in this sample using three different models: a sin- 



gle Sersic, a deVExp and a SerExp (see Meert et al. 2012a,b 
for details). We then use the best-fit parameters from these 
different fits to generate three synthetic images for each ob- 
ject. In this way, we have, in effect, three different mock 
SDSS catalogs (see Meert et al. 2012a for detailed tests). If 
galaxies were, in reality, e.g. two-component deVExp mod- 
els, then only our deVExp mock catalog would be realistic 
- performing profile fits (e.g., using the other two models) 
to this catalog should return results which are similar to 
those when fitting to the SDSS data. Moreover, although 
all three catalogs will contain correlations between n, total 
luminosity, half-light radius, etc., these correlations are only 
guaranteed to be like those in the SDSS data for this (in this 
case, deVExp) mock catalog. 
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Figure B3. Parameters n ser of the bulge and B/T obtained from 
fitting the two-component SerExp model to mock galaxies gener- 
ated using input deVExp model (left) and SerExp model (right). In 
the panel on the left, the fits correctly return values of n sor ~ 4; 
in the panel on the right, the distribution resembles the input 
one: notice that this one indicates that bulges do not necessarily 
have n = 4. 



B2. 1 Fitting to a profile which is truly a single Sersic 

We begin with the case in which PyMorph is asked to fit what 
is in reality a single Sersic profile of index n (i.e. we use the 
mock galaxies generated using a single Sersic profile) with a 
single component Sersic, and with deVExp and SerExp pro- 
files. The distribution of input n values used to simulate the 
mock galaxies is that which one obtains from fitting single 
Sersics to the parent (magnitude limited) sample. Rather 
than showing the fits themselves, we present our results in 
the parameter space of the best-fit n versus best-fit B/T. In 
all cases, darker shading indicates regions in the parameter 
space that are more heavily populated. 

The top row in Figure IBT1 shows results for input single 
Sersic mock galaxies. The two panels on the left show B/T 
values determined from the deVExp fits, and the two on the 
right are from SerExp fits. We describe the deVExp results 
first. The top left panel of Figure IBTI shows the distribution 
of the sample in best fit n— B/T space, and the next panel 
to the right shows the result of restricting the analysis to 
narrow ranges of input n. The different colors show the dis- 
tribution in fitted n and B/T for input n in the range — 2, 
2 — 4, 4 — 6 and 6 — 8 (we show the regions which enclose 
25%, 50% and 75% of the points). Comparison with the val- 
ues along the y-axis shows that PyMorph correctly returns 
the input n values. 

The distribution in the n— B/T plane is clearly non- 
trivial. For n < 4 there is a tight correlation between the 
value of n returned by the single component and B/T from 
the deVaucouleurs-exponential fit: B/T— > 1 as n — > 4. But 
as n increases beyond 4, B/T begins to decrease again. I.e., 
B/T is not a monotonic function of n. Since the deVExp 
profile only has n — 1 or n = 4 components, to fit n > 4 
profiles PyMorph requires more and more of an exponential- 
like component, i.e. B/T decreases. (The figure does not 
show this, but the fit returns bulge half-light radii which 
are ever smaller fractions of the half-light radius of that of 
the input Sersic profile.) As a result, for 1/2 <B/T< 1, 
the distribution of n at fixed B/T appears bimodal. This 
shows that, unless one is certain that large values of n do 
not occur in nature, then, especially around B/T~ 0.7, B/T 
values may be misleading, if not meaningless. 

The two panels on the right show the corresponding dis- 
tribution for SerExp; they are clearly different from those for 
deVExp. This is primarily because PyMorph correctly assigns 
the entire profile to the bulge (Sersic) component, except 



when the input n ~ 1, since then which of the two n = 1 
components should be called the bulge is ambiguous. (We 
have checked that, when n ~ 1 and B/T < 1, then the half- 
light radius of the 'bulge' component is indeed the same as 
that of the total: i.e., the two components differ only by the 
value of B/T.) The fact that B/T is not exactly equal to 
unity is a measure of the error in B/T which comes from 
the extra degree of freedom associated with having a second 
component with which to fit the profile. 

B2.2 Fitting to a profile which is truly a deVExp 

The second row shows results when the input profile used 
to simulate the mock galaxies is a two component deVExp 
model (the distribution of input B/T values is obtained from 
fitting deVExp models to the SDSS parent magnitude limited 
sample) . This two-component profile is then fit with a single 
Sersic to get n; B/T comes from fitting a deVExp model (two 
panels on left) or a SerExp model (two panels on right). The 
overall (grey-scale) distributions are rather different than in 
the corresponding panels in the top row. This is the first 
hint that the distribution of fitted n-B/T can be used as 
a diagnostic of the true profile shape. Different colors show 
results for narrow bins in input B/T; these indicate that 
PyMorph indeed returns the correct values when it fits the 
right model. The additional freedom when fitting a SerExp 
profile to what is really a deVExp means that, in the panel 
on the far right, the distribution of fitted B/T at fixed input 
B/T is slightly broader than when fitting a deVExp. 

B2.3 Fitting to a profile which is truly a SerExp 

Finally, the bottom row shows results when the input model 
used to simulate the mock galaxies was a SerExp (with n and 
B/T values chosen from fitting the SDSS parent sample to 
a SerExp model). The results here differ from those in the 
row above in subtle ways, perhaps most appreciably in the 
upper right corner (large fitted n and B/T) of the bottom 
right plots. 

In this case, we also show (Figure [B3} the ntuige— B/T 
plane, where both ribuige and B/T come from fitting a 
SerExp model to mock images generated using input deVExp 
(left) and SerExp (right) profiles. The panel on the left shows 
that PyMorph correctly returns n ~ 4 when it should; we 
have checked that the distribution in the panel on the right 
is similar to the input one, again suggesting that PyMorph is 
working well (Meert et al. 2012a). 

B2.4 Fitting to SDSS images 

Figure IB2I shows a similar analyis of SDSS images. In the 
two panels on the left, n comes from fitting a single compo- 
nent Sersic, and B/T from fitting a two-component deVExp. 
In the panels on the right, B/T comes from fitting a two- 
component SerExp. Notice that the gray scale plots are very 
unlike those in the top row of Figure [Bll and most like those 
in the bottom row. This suggests that SDSS galaxies are al- 
most certainly not single-component systems. 

In addition, of the two-component models, the SerExp 
model appears to be more like the data than is the deVExp. 
This is because, when B/T comes from fitting a SerExp, 
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Figure B4. Fitted n ser vs fitted B/T for simulated galaxies that 
are assumed to be single Sersic profiles (top), two-component 
deVExp profiles (middle) and two-component SerExp profiles (bot- 
tom). In all cases, the y-axis shows n scr returned by fitting a 
single Sersic profile to the image. In the left column, B/T is 
obtained from fitting a two-component deVExp model; the right 
column, B/T is determined from fitting a SerExp model. The 
density is shown in four bins colored by output absolute mag- 
nitude: -24 < M r < -23 (red), -23 < M r < -22 (green), 
-22 < M r < -21 (cyan), -21 < M r < -20 (blue). 
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Figure B6. Comparison of total apparent magnitude (top) and 
luminosity (bottom) returned from single Sersic and SerExp fits 
to simulated SerExp (left) and real SDSS (right) galaxies. The 
error bars show the la rms scatter around the median. 
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Figure B5. Similar to Fig. IB4l but for real SDSS galaxies 



then the SDSS data (third panel from left) populate the 
large n— B/T corner which input SerExp models also fill, 
but input deVExp models do not (c.f. Figure iBl) . There is a 
more subtle difference when B/T comes from the deVExp fit 
(left most panels) in Figures IB 1 1 and IB2I the SDSS shows 
a rather well-defined ridge at the boundary of the large 
n— B/T corner, which appears to be more separated from 
the peak at small n. This separation is more apparent for 
the input SerExp models than for input deVExp. 

Since we cannot classify the objects by the true value of 
n or B/T, the colors (contours in Fig. IB2|I show the result 
of restricting the analysis to objects which are most likely 
to be ellipticals (red) to least likely (blue), as determined by 
Huertas-Company et al. (2011). This shows that the ellipti- 



Figure B7. Comparison of angular (top) and physical (bottom) 
half-light radii returned from single Sersic and SerExp fits to sim- 
ulated SerExp (left) and real SDSS (right) galaxies. The error 
bars show the lcr rms scatter around the median. 



cals do indeed have large values of n, and spirals the lowest, 
as expected. 

To provide a slightly more straightforward compari- 
son between simulations and data, we have considered the 
n — B/T distribution for objects in narrow bins in (output) 
luminosity. Figures IB4I and IB5I show results in simulations 
(the same fits used for Figure IB 1 jl and in the SDSS (cf. 
Figure IB2[I , respectively. These too indicate that the two- 
component models are more like the data, with the SerExp 
marginally favoured (the two panels in Figure lB5l look more 
like the bottom than the middle panels of Figure IB4[) . 



B3 Biases from fitting single Sersic profiles 

The analysis above shows that a single component Sersic 
profile is not as good a description of SDSS galaxies as one 
with two-components. Since such single component fits are 
much simpler to perform, and are commonly used, it is inter- 
esting to ask if they lead to significant biases in commonly 
used parameters. E.g., one might expect the total light to 
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Figure B10. Biases in the estimated luminosities and sizes which come from fitting single Sersic and two-component SerExp profiles to 
images which are really pure Sersics (left) and two-component SerExps (right). 
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Figure Bll. Biases in the estimated luminosities (left) and sizes (right) of the total (green), bulge (magenta) and disk (blue) components 
in SerExp fits to SerExp images. The estimated total and disk components are usually unbiased, whereas the bulges tend to be too big 
and too bright at large luminosities. 
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Figure B8. Comparison of luminosities returned from single Ser- 
sic and SerExp fits to real SDSS galaxies, color coded by best-fit 
n (top), best-fit B/T (middle), and p(Ell) (bottom). 



Figure B9. Comparison of angular half-light radii returned from 
single Sersic and SerExp fits to real SDSS galaxies, color coded 
by best-fit n (top), best-fit B/T (middle), and p(EU) (bottom). 



be a reasonably robust quantity, so different models for the 
shape of the profile may still return consistent values of Ltot- 
The top left panel of Figure IB6I shows that the apparent 
magnitudes returned by single Sersic fits to the objects in 
our SerExp mock catalog are quite accurate, with a ten- 
dency for the Sersic fits to return an overestimate by about 
ten percent at the bright end. The top right panel shows 
that a similar comparison for the objects in the SDSS pro- 
duces similar results. The bottom panels show the impact 
of these small biases on the inferred luminosities. 

Figure [B7l shows a similar analysis of the half-light radii: 
the single Sersic fit tends to overestimate the sizes by about 
ten/fifteen percent, particularly for the largest objects. The 



largest and/or most luminous galaxies tend to have large 
n and/or intermediate to large B/T. Therefore, this bias is 
worst for objects that are likely to be ellipticals. We show 
this explicitly in Figures [B8I and TB9I 

We have presented a novel diagnostic of whether or not 
the surface brightness profiles of galaxies are better thought 
of as having one or two components. The method works 
by fitting a number of single and two-component models to 
the image, and then studying the distribution in the n-B/T 
plane defined by the Sersic index n associated with the sin- 
gle component fit and the ratio B/T of bulge to total light 
in the two component fit. The way SDSS galaxies populate 
this plane suggests that they are not single component Sersic 



© 0000 RAS, MNRAS 000, 000-000 



22 Bernardi et al. 



systems. Rather, their distribution in rt-B/T is more similar 
to that expected of two-component systems, with a Sersic 
+ exponential model faring somewhat better than the tra- 
ditional deVaucouleurs bulge with exponential disk model 
(Figures IB II and IB2[) . I.e., in bulge dominated systems, al- 
lowing n / 4 provides a significantly improved fit. Indeed, 
we even find bulges with n > 4 in the SDSS (Figure |B3[) . 

Our conclusion that the SerExp model is preferred is 
consistent with a recent analysis of the MGC, indicating 
that at least half of the galaxies at z ~ 0.1 are two compo- 
nent SerExp systems (Allen et al. 2006). Forcing the light 
profiles to be fit by a single Sersic profile leads to biases in 
the inferred total luminosities and half-light radii of galax- 
ies: both are overestimated by about ten percent, especially 
at large luminosities and sizes. The bias is dominated by 
objects with large fitted n and/or intermediate values of fit- 
ted B/T; these are objects that are likely to be early-type, 
but have a significant exponential component, so the as- 
sumption of a single profile is particularly bad. In contrast, 
objects that are likely to be late-type are unbiased. 

These biases have a small systematic effects on the size- 
luminosity correlation of objects that are likely to be early- 
types, and are presented in the main text. 



APPENDIX C: CORRELATION BETWEEN 

-Rbulge/ Rh OR Roisk/Rh AND B/T 

Figure [TT] of the main text showed a correlation between the 
bulge(disk) to total size and B/T. This correlation arises 
because the Sersic profile is 



Z(r) = /„ exp[-(r/r„) l/B 



(CI) 



so the ratio of the light within r to the total light in the 
profile is 



L n (< r) _ JZ /rn dxx exp(-a 1 /") 



J °° dxx exp(— x 1 



72n [0,(r/r n ) 1 ''' 1 ] 



(C2) 

where 72„ is the incomplete Gamma function. (For integer 
n, it can be written in terms of exp[— (r/r„) llln ] times a 
polynomial of degree In - 1 in (r/r n ) 1/n .) Therefore, the 
half-light radius of a SerExp profile satisfies 



B 



' 72n 



0, 



1/n 



+ 1 



7-> 



0, 



l/n 



(C3) 

For a given B/T, the right hand side is a function of rh/r n 
and Th/ri = (rh./r n )(r n /ri), so it defines a different curve 
for each r n /ri, where r n — rb u i gc /(1.992n — 0.327) and 
n = r<iisk/l-67. Note that the curves are independent of 
luminosity L; therefore L dependence only enters if the dis- 
tribution of r n /ri and/or B/T depend on L. 

Figure [Cll shows -Rbui g e/-Rhi as a function of B/T for 
the early- (top) and late-type (bottom) samples defined in 
the main text for galaxies with —21.5 > M r > —22.5; re- 
sults are similar at other luminosities. The curves show the 
predicted relations (equation IC3|) for a deVaucouleur bulge 
(n = 4) with exponential disk. These depend on the ratio 
fl!disk/-Rbuige, for which we have chosen 2, 4 and 6. 

The top panel shows a very strong correlation between 
flW g o/-Rhi and B/T (at this fixed M r ) for the early-type 
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Figure CI. Correlation between i?bulgc/-f?hl and B/T for early- 
types (top panel) and between -Rdisk/^hl an d B/T for late-types 
(bottom panel). Although only objects with —21.5 > M r > —22.5 
are shown, we see qualitatively similar behaviour at other lumi- 
nosities. Dashed, dotted and solid curves show the expected cor- 
relation for n = 4 bulges with Rdisk/^bulgc = 2, 4 and 6. 



sample. Clearly, if 20% of the light is in a disk component, 
then the size is affected by at least this fraction. The well- 
known correlation between L and B/T, and the fact that 
early-types span a large range of B/T, means that the bulge 
and early-type size-luminosity relations can be quite differ- 
ent. It is perhaps surprising that the half-light radius of the 
disk component is typically more than 3-5 times larger than 
that of the bulge, particularly at B/T< 0.7. We return to 
this shortly. 

The bottom panel shows i?disk/-Rhi and B/T for 
the late- type sample. Most of the sample has B/T< 
0.2. Comparison with the smooth curves indicates that 
fldisk/^buigc ~ 5 for most of the sample. In this case we do 
expect the disks to be substantially larger than the bulges, 
so the results are sensible. 

Figure [TT] in the main text shows this same information 
in a different format, which allows for a more direct under- 
standing of the impact this correlation has on the relations 
shown in Figure [5] And Figure in the main text shows 
that -Rdisk/^buigc is indeed substantially larger than unity. 

To address the question of large Rd/Rb in our early- 
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type sample, particularly at smaller B/T, Figures [C2l and [C3l 
show two examples. The format in both cases is the same: 
The top left panel shows a ~ 20 arcsec field centered on the 
object, to get an idea of whether or not the object is in a 
crowded field. The top right panel provides a closer look at 
the object. The panel just below it shows the best-fit Ser- 
Exp model, and the middle left panel shows residuals from 
this fit. The bottom left panel shows the one-dimensional 
surface brightness profile, and our Sersic (solid magenta), 
deVExp (solid blue) and SerExp (solid red) fits; dotted and 
dashed curves show the corresponding disk and bulge com- 
ponents. Bottom right panel shows the associated residuals. 
The legend along the left shows the values of many quanti- 
ties returned by the fits, and other information, such as the 
BAC p(type), for the object. 

The object in Figure IC2l B/T= 0.71 and R d /R b ~ 10 
is very likely to be an elliptical: p(Ell)= 0.87. The Sersic 
and SerExp fits return almost the same magnitudes (M r ~ 
—22.2) and total half light radii (~ 3.15"). However, n = 
7.15 for the single Sersic fit, but n — 4.79 for the SerExp 
bulge. For the SerExp, as for the deVExp fits, the second 
component is clearly necessary. The Xdof values for these 
fits are similar. In Figure [C3l the Sersic magnitudes and half 
light radii are slightly larger, but otherwise the qualitative 
trends are the same: the single Sersic fit requires large n, 
and the second component in the SerExp fit clearly requires 
large R d /R b ~ 6. 
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Note: The Id data is calculated using background-subtracted data. The 2d data is shown with background included. 



Figure C2. Example of an early-type galaxy with —21.5 < M r < —22.5, large R^/Rb and B/T~ 0.7. Top left panel shows a ~ 20 arcsec 
field centered on the object; top right panel provides a closer look. Middle right panel shows the best-fit SerExp model; middle left panel 
shows residuals from this fit. Bottom left panel shows the one-dimensional surface brightness profile (symbols), and our Sersic (solid 
magenta), deVExp (solid blue) and SerExp (solid red) fits; dotted and dashed curves show the corresponding disk and bulge components. 
Bottom right panel shows the associated residuals. Legend along the left shows the values of many quantities returned by the fits, and 
other information, such as the BAC p(type), for the object. 
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Note: The Id data is calculated using background-subtracted data. The 2d data is shown with background included. 



Figure C3. Same as previous figure, but for another early-type galaxy selected at random from among those with the same M r and 
B/T range. 
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